Width of exotics from QCD sum rules : tetraquarks or molecules? 
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We investigate the widths of the recently observed charmonium like resonances X(3872), Z(4430) 
and ^2(4250) using QCD sum rules. Extending previous analyses regarding these states as diquark- 
antiquark states or molecules of D mesons, we introduce the Breit-Wigner function in the pole term. 
We find that introducing the width increases the mass at small Borel window region. Using the 
operator-product expansion up to dimension 8, we find that the sum rules based on interpolating 
current with molecular components give a stable Borel curve from which both the masses and 
widths of these resonances can be well obtained. Thus the QCD sum rule approach strongly favors 
the molecular description of these states. 

PACS numbers: 11.55Hx,12.38.Lg, 12.39.-x 



I. INTRODUCTION 

Since the discovery of a charmonium like resonance 
X{3872) by Belle Collaboration plenty of similar res- 
onances have been observed in the decay of B mesons. 
In particular, the recently observed Z+(4430) in ir^ip' in- 
variant mass spectrum has a charge 2] , and consequently 
can not be a simple charmonium. There arc already a 
number of interpretations on the structure of these res- 
onances 0. Although the nature of these states is still 
an open question, tetraquark state and molecular state 
are both intriguing possibilities. Since the decay prod- 
ucts are a charmonium and a pion, it is natural to expect 
that the parent contains four quarks including c and c. 
On the other hand, the masses of X (3872), Z+(4430) and 
the recently observed resonancelike structure ^^(4250) 
are very close to thresholds of two Z?-meson states, 
D*D, D*Di and DiD, respectively, and, therefore, it 
is very tempting to interpret these states as molecular 
states. 

Motivated by these facts, QCD sum rules (QCDSR) 
have been extensively used to study these resonances. In 
Ref. 5'], X{3872) was analyzed by assuming it to be a 
J-^*-^ = 1++ tetraquark (ccqq) state. Although the re- 
sult agreed with the experimental data, an analysis us- 
ing a current composed of a D*D molecule shows bet- 
ter operator-product expansion (OPE) convergence and 
closer agreement with experimentally observed mass 
strongly suggesting a molecular nature of X(3872). 

In Ref. d, Z+(4430) was considered as a D*Di 
molecule and a good agreement with data was obtained, 
while the tetraquark description has been found to be 
unsatisfactory [8']. Similarly, we have applied the molec- 
ular description to the most recent data of Z]^(4050) and 
(4:250) 0], and found that a DiD molecular state can 
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be attributed to Z^(4250). However, it was not possible 
to explain Zj^(4050) as a molecular D*D* state Q. In 
all of the calculations above, however, small but finite 
width was not taken into account, in spite of that fact 
that Z+(4430) and ^2^(4250) have widths r^+(443o) = 
45t\ltll MeV i and T^^^,,,,^ ^^ I77tlttll' MeV i, 
respectively. Of course these widths are much smaller 
than their masses, around 4 GeV. However, the effect 
of the width should be examined in order to clarify the 
structure of these states. 

In this paper, we extend the previous analyses @, 
0, H, [§] to include the effect of finite width and give fur- 
ther consideration on the possibility that these states can 
be considered as tetraquarks or molecules. The width is 
usually not calculated in QCD sum rule approaches as 
the OPE are usually restricted to three terms; perturba- 
tive, dimension four and dimension six terms. Hence the 
phenomenological sides are composed of three unknown 
parameters; mass, continuum, and overlap constant. In 
the present analysis, the OPE are composed of operators 
with four different dimensions. Therefore, an analysis 
including the width is sensible. 

In the next section, we give a brief review of our 
QCDSR analyses. In Sec. IIIIl we discuss some general 
features of effect of finite width. Quantitative analyses 
of the exotic states are given in Sec. IIVI Section |V] is 
devoted to the summary. 



II. QCD SUM RULES 

The QCD sum rules for mesons are based on the two- 
point function of a current j{x) describing a desired state 

n(g)=z I dW?-(0|Tb-(x)jt(0)]|0), (1) 

and the dispersion relation 

11(9) / ds -\- (subtraction terms) (2) 
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with p{s) being the spectral density. While computing 
the two-point function in terms of quarks and gluons by 
making use of operator product expansion (OPE), which 
takes into account non-perturbative effect through QCD 
condensates, one models the hadronic spectral density 
^phcn^g^ with a pole describing the ground state and a 
continuum, namely. 
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In the narrow width approximation, the pole part of 
the hadronic spectral density is set to a delta function 
pP°^'^{s) — X^S{s — m?), with (O|j|meson) = A being the 
overlap of the current and the physical meson. In this 
work, we replace this part with the relativistic Breit- 
Wigncr function to take the width into account. The 
continuum part above the threshold so is given by the 
result obtained with the OPE 



p™"*(s)=pOP^=(,sM,s-so), 



(4) 



with 0{x) being the step function. The OPE side is cal- 
culated up to leading order in as and condensates up to 
dimension eight are considered. Currents and OPE terms 
used in this work are taken from Ref. [E B S H Q • The 
correlation function in the OPE side can be expressed as 
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where p°^^(s) = 7r~^Imn'-*^^(s). Then, we can extract 
the pole term by equating the OPE expression and the 
phenomenological expression, making the Borel transfor- 
mation on both sides and then transferring the contin- 
uum contribution ^ to the OPE side. The sum rule is 
hence given by 



-Hn™<9«>(M2). (6) 



Note that the left-hand side becomes A^e"™ in the 
narrow width approximation. In this case, one can ex- 
tract the pole mass by taking the ratio between the 
derivative of Eq. ^ with respect to 1/M^ and Eq. Q 
itself. In the present work, we introduce the width by 
employing the Brett- Wigner function to the pole contri- 
bution 
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The mass and width are determined by looking at the 
stability of mass against varying Borel mass M^, as usual. 
The relevant Borel window is determined by the conver- 
gence of the OPE for the minimum M^^^ and the pole 
dominance criterion for the maximum M^^^. As usual, 
we determine M^^^ by requiring the dimension eight con- 
densate contribution to be less than 10% of the total 



OPE and M^^^ from more than 50% pole dominance. 
We calculate the mass by fixing a width and solving the 
equation for the ratio 
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where Ii'^^^{M'^) is the right-hand side of Eq. 
pP°i'=(s) is given by Eq. ^. 
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FIG. 1: (color online) . Right-hand side of Eq. (|8]) as a function 
of Borel mass with fixed m and F. Upper and lower 
panels stand for the case of m = 4 GeV and m = 4.4 GeV, 
respectively. 



III. GENERAL FEATURES 

Since the left-hand side in Eq. ([5]) does not change by 
introducing a width, one can investigate how the mass 
changes only by looking at the behavior of the right-hand 
side. 

Figure [T] shows the right-hand side of Eq. ([8]) as a func- 
tion of Borel mass IvP . One can see that it is a monotonic 
function of if both m and F are fixed. It rapidly in- 
creases at small and then asymptotically reaches to 
the Breit-Wigner mass m. In analyses of QCDSR, we 
solve Eq. ^ with a value of left-hand side of Eq. ([5]) 
given by the OPE side and the continuum. This pro- 
cedure corresponds to finding an intersection between a 
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the top panel. This means that if one gets the mass 4200 
MeV in the T = case with stabihty at = 2.5 GeV^ 
the maximum width of this state is hmited to 50 MeV. 
This also gives the constraint on possible mass and width 
values. 
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FIG. 3: (color online.) Results for X(3872) as a tetraquark 
state with width. Continuum thresholds so is taken to be 
^5^= 4.2 GeV. 



FIG. 2: (color online.) Same as Fig. (TJ but as a function of 
the Breit-Wigner mass and for a fixed Borel mass. 



horizontal line denoting the value of mass in the F = 
and the curves of a given Borel mass in the figures. For 
example, if one has to = 3.8 GeV in the F = case, a 
possible solution at = 2.0 GeV^ and m = 4 GeV is 
F ~ 40 MeV. If one sets m = 4.4 GeV, F ~ 60 MeV is one 
of the possible solutions. The best solution is determined 
by looking at the stability against M^. From the mono- 
tonic behavior seen in Fig. [Jl one notes that introducing 
the width increases the mass especially at small Borel 
mass region. Hence, if one gets the mass in the F = 
case which monotonically increases as increases, it 
will be improved by including the width. This fact gives 
a guideline on the QCDSR analyses.^ 

We also plot the direct relation between the Breit- 
Wigner mass and the mass in the F = case in Fig. [51 
Here we fixed the Borel mass = 2.5 GeV^ in the top 
panel and 3.0 GeV^ in the bottom panel, which are typi- 
cal values satisfying the stability criterion in the QGDSR 
analyses below. One can see that deviation from the mass 
in the F = is larger for larger mass to and smaller Borel 
mass M^. One should note that it is no longer monotonic 
as a function of to at smaller and large F, as seen in 



^ This is not a completely general result of QCDSR; as seen in 
[ld |. introducing width loads to smaller mass when is large 
compared to m. 
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FIG. 4: (color online). Results for X(3872) as a D* D 
molecule. The crosses indicate lower and upper limit of the 
Borel window, respectively. 



IV. RESULTS 

For parameters in the QCDSR analyses, we use the 
same parameter set as in the previous works and assume 
the factorization of the higher dimensional condensates. 
Namely, = 1.23 GeV, (gg) = -(0.23GeV)3, (£,2^2) = 
0.88GeV^ {qga-Gq) = ml{qq) and = O.SGeV^ Here 
we ignored the possible uncertainties in these parame- 
ters. Because the uncertainties are those related to the 
OPE side, which are unchanged by including the width. 
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possible errors on masses will not be so different from the 
ones estimated in the previous works. 

From the consideration in the previous section, one 
finds that tetraquark configurations for Z(4430) exam- 
ined in Ref. P are ruled out even if the width is taken 
into account. In the = 0^ results (Fig. 3 of Ref. [1]), 
the mass in the F = case shows a good stability. In- 
corporating width does not improve the stability, but it 
raises the value of mass, which is already bigger than the 
experimental value in the F = case. If one assumes 
= 1~ tetraquark state, it is shown that the mass is 
more than 300 MeV larger than the experimental value, 
and the functional behavior with respect to is mono- 
tonically decreasing Q . Both of these features are only 
worsened by introducing the width in the calculations. 
The failure to explain the Z(4430) with tetraquark struc- 
tures cannot be corrected by taking width into account. 
In the case of X{3872), the value of mass in the F = 
limit is found to be in agreement with the data Q while 
functional behavior against is monotonically decreas- 
ing. In this case, however, the experimental width has 
been found to be around 2.3 MeV Then, such a 
small width hardly improves the Borel curve, as shown 
in Fig. W 

Contrary to the sum rules obtained with interpolating 
currents based on large tetraquark component, molecular 
descriptions are much more promising. First, let us con- 
sider a counter example of Fig. [3l The current for D*D 
molecule is obtained by replacing the strange quark by a 
light quark in Ref. 

where a and b are color indices. The combination 
- D^D*^ has J^<^ = 1++ as the X(3872) meson. 
To show that this current has positive charge conjugation 
we notice that under charge conjugation transformation 
the two terms in Eq. ^ transforms as: 

C[iUal5Ca)icblf,Ub)]C^'^ = - [(CaTS^a) (ufcT^Cb)] 
C[{Ca"f5Ua)iub"ftJ.Cb)]C^^ = -[{Ua"f5Ca){cbJfj.Ub)]ilO) 

Therefore, one obtains 

Cj^C-'^j^. (11) 

The symmetrical combination, D'^D*'^ + D^D*^ , would 
provide exactly the same mass, within our sum rule ap- 
proach. 

The resultant OPE expressions for D*D molecule are 
obtained by putting = and (ss) = {qq) in the ex- 
pressions in Ref. [6,] . Figure S] shows the result of mass 



^ In Ref. 0, the mass of X(3872) was evaluated by including con- 
densates up to dimension five and higher dimensional ones were 
used to estimate errors. In this paper, we have included all the 
condensates given in Ref.[||. 



of the D*D molecular state with F = 0,2,4,6,8 MeV. 
One sees that the molecular description gives better sta- 
bility than the tetraquark case. Furthermore, although 
effect of width is not larger, we can fit the experimental 
mass 3872 MeV and width F < 2.3 MeV simultaneously 
with a continuum threshold ^/sq — 4.38 GeV, and ob- 
tain a broad Borel window. Consequently, it is strongly 
favored that X(3872) is a D*D molecular state. 

Next, we consider Z(4430) states as a D*Di molecule. 
The current and OPE expressions are given in Ref. 0]. 
It has been shown that the molecular description gives a 
mass which agrees with the experiment well. Since the 
mass in the F = case is monotonically increasing func- 
tion of M^, it is expected that incorporating the width 
improves the stability according to the result shown in 
SecHIl 

Figure [5] shows the result of D*Di molecule with var- 
ious continuum thresholds. In Ref. the continuum 
threshold is determined as = 4.8 — 5.0 GeV. The 
center value = 4.9 GeV case is plotted in the right- 
bottom panel. In this case, the mass in the F = case 
agrees well with the experimental value. As introducing 
the width raises the mass, however, the mass becomes 
larger than experiment when F is as large as the experi- 
ment. One notes that the stability becomes much better 
when F ~ 30 MeV. Other three panels show the cases 
with lower continuum thresholds. Especially ^^sq = 4.6 
GeV case reproduces both mass and width quite well. 
One can see that F ~ 40 MeV gives the best stability of 
the mass, which perfectly agrees with the the experiment. 
In the lower continuum thresholds case, however, we have 
to relax the criterion for the allowed region of sum rule 
analyses, i.e., Borel window. Since the mass of Z(4430) 
is close to D*Di threshold, it might be plausible that 
the continuum contribution becomes larger when D*Di 
forms a molecule, as assumed in this calculation. The ar- 
rows in the figure indicate the Borel masses determined 
from various values of the relative contributions of the 
dimension eight condensates (for M^j^^, upward arrows) 
and the continuum contribution (for M^^^, downwards 
arrows). One sees that reasonable Borel windows open 
if one relaxes the condition for either M^- or M^„^, or 

mill iiicix ' 

both of them. 

One notes that there is a truncation of the curves in 
each panel, especially for large width data. This is due to 
the nature of the Breit-Wigner function, shown in upper 
panel of Fig. [21 that the right-hand side of Eq. ([S]) has 
the maximum at low Borel mass and large width region. 
This fact appears as an absence of the solution of Eq. [H] 
for a fixed Borel mass and F. Hence, it expresses a maxi- 
mum width allowed by the QCDSR for each value of the 
continuum threshold. 

Finally we consider the recently discovered Z^(4250) 
a,s a DiD molecule. The current and resultant OPE ex- 
pressions are given in Ref. In Ref. 0, it is shown 
that DiD molecular description gives a reasonble agree- 
ment with the ^2^(4250) mass. As mentioned above, this 
state has a large width whose effect should be examined. 
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FIG. 5: (color online). Results for D* Di molecule state. Each panel shows different continuum threshold case. Upward and 
downward arrows indicate the region of the Borel window M^j^ and M^^^, respectively. Associated numbers in percent denote 
the dimension eight condensate contribution for upward arrows and continuum contribution for downward ones. 



Indeed, dependence of the mass looks promising be- 
cause it is a monotonically increasing function of as 
in the case of D*Di molecule. 



Figure [H] shows the results for the DiD molecule. As 
in Fig. [51 we plot the mass for various widths in each 
panel. The continuum threshold = 4.6 GeV corre- 
sponds to the center value in the previous analysis Q. 
One can see that the stability is achieved in all the cases. 
In the largest ^/so case, the optimized width is around 
40 MeV, which is a little smaller than the experimental 
value. Reducing the continuum threshold leads to larger 
width; for ^ = 4.5 GeV, we obtain the F = 60 MeV. 
Up to this value, no relaxation of the Borel window crite- 
rion is needed. Taking ^/sq < 4.4 GeV makes the width 
much closer to the experiment, however, we need to relax 
the condition for the dimension eight condensates and/or 
the continuum contribution as in the case of Z(4430) to 
validate QCDSR. In such cases, we can obtain F ~ 80 
MeV (y/s^ = 4.4 GeV) and F ~ 100 MeV = 4.3 

GeV). Hence, our analyses support an existence of DiD 
molecular state with large width and its manifestation as 
2'2'"(4250) resonance. 



V. SUMMARY 

In summary, we have extended previous QCDSR anal- 
yses of exotics to include the total width, by employing 
the Breit-Wigner function to the pole term. As a gen- 
eral feature, for the cases where the predicted mass for 
F = increases with increasing Borel mass M^, introduc- 
ing the width increases the predicted mass at small Borel 
mass region, and improves the Borel stability, From this 
point of view, none of the sum rules based on interpolat- 
ing currents with tetraquark components are favored. On 
the other hand, sum rules based on interpolating currents 
with molecular description as D*D, D*Di and DiD, are 
shown to give valid sum rules for X{3872), Z+(4430) and 
Z^{4250) respectively. For X(3872), the inclusion of the 
width slightly modify the mass, leading to a better agree- 
ment with the experimental result. For Z+(4430) and 
^^(4250), molecular description proposed in Refs. 0,0 
are largely improved by introducing the width. We have 
obtained stable results with Tz{4i3o) — 40 MeV and 
rz2(4250) — 40 — 100 MeV. These results strongly sup- 
port the previous results based on F ^ limit, that 
the Z+(4430) and .Z^(4250) resonances are strong can- 
didates for molecular states. Moreover, we have estab- 
lished that the QCD sum rule with four OPE terms has 
sufhciently rich structure so that an estimate of the total 
width is also possible. 
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FIG. 6: (color online). Results for DiD molecule. Symbols are similar to Fig. [5] 
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